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Abstract: Coannihilation processes provide an important additional mechanism for re- 
ducing the density of stable relics in the Universe. In the case of the stable lightest neu- 



tralino of the MSSM, and in particular the Constrained MSSM (CMSSM), the coannihi- 
lation with sleptons plays a major role in opening up otherwise cosmologically excluded 
ranges of supersymmetric parameters. In this paper, we derive a full set of exact, ana- 
lytic expressions for the coannihilation of the lightest neutralino with the sleptons into all 
two-body tree-level final states in the framework of minimal super symmetry. We make no 
simplifying assumptions about the neutralino nor about sfermion masses and mixings other 
than the absence of explicit CP-violating terms and inter-family mixings. The expressions 
should be particularly useful in computing the neutralino WIMP relic abundance with- 
out the approximation of partial wave expansion. We illustrate the effect of our analytic 
results with numerical examples and demonstrate a sizeable difference with approximate 
expressions available in the literature. 
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1. Introduction 



The relic density of stable weakly-interacting massive particles (WIMPs) is determined 
primarily by how efficiently their number density in the early Universe can be reduced. 
In the case of the most popular WIMP candidate: the lightest neutralino of minimal 
supersymmetry (SUSY), assumed to be the lightest SUSY particle (LSP), there are two 
generic mechanisms |2|. First, the neutralino can pair-annihilate into ordinary particles. 
Second, in some cases they can coannihilate with some some other species if these are 
nearly mass-degenerate with the LSPs. 

The standard mechanism of neutralino pair-annihilation has been considered in much 
detail in many papers 0. In particular, complete sets of neutralino annihilation cross 
sections were provided in ||] (see also Q) in the case of partial wave approximation. Exact 
expressions applicable both near resonances and new final-state thresholds were recently 
published in Q. 

The mechanism of coannihilation was originally pointed out by Griest and Seckel p|. 
It applies when there exists some other species x' which is not much heavier than the 
WIMP and may therefore be still present in the thermal plasma at the time of WIMP 
decoupling. Coannihilation becomes important if its annihilation with the WIMP (and/or 
itself) is equally, or more, efficient than the pair-annihilation of the stable WIMPs. 

These circumstances in particular are naturally realized in the case of the higgsino-like 
lightest neutralino LSP in the framework of the Minimal Supersymmetric Standard Model 
(MSSM). In this case the next-to-lightest neutralino and the lightest chargino are almost 
mass-degenerate with the LSP || [7j. In fact, the coannihilation in this case is so efficient 
that it has a devastating effect on the relic density of the higgsino-like neutralino, which 



is the type of LSP strongly disfavored by naturalness |8[ and mass-unification y, |10|, |ll[. 
It even plays some role in the strongly prefered case of the bino-like LSP ||. 

Since, in the framework of general softly-broken low-energy SUSY, scalar superpartner 
masses are a priori unrelated (or at best loosely related) to the neutralino sector, in principle 
one might assume that any of the scalar superpartners could be light enough to participate 
in coannihilation with the neutralino LSP - the case that would be technically rather 
challenging and in any case not particularly well-motivated. 

A more realistic scenario is the one in which one of the scalar partners of the top- 
quark or the r-lepton is rather light and nearly degenerate in mass with the neutralino 
LSP. This is because the off-diagonal elements in their 2x2 mass matrices can under some 
circumstances greatly reduce one of the eigenmasses relative to the other. The case of 
neutralino-stop coannihilation was considered in [12] in the framework of the MSSM and 
recently shown in [13] to be also applicable in the case of the Constrained MSSM (CMSSM) 
but only for rather large values of the trilinear soft SUSY-breaking term Aq. 

The importance of the neutralino coannihilation with the lighter of the two staus in 
the framework of the CMSSM was pointed out in [14]. In this model, when the common 
gaugino mass parameter my 2 is much larger than the common scalar mass parameter mo, 
it is the lightest stau that is the LSP [11]. The effect of the neutralino-stau coannihilation 
is to open up a narrow corridor [14] just above the boundary of equal neutralino-stau 
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masses into an otherwise cosmologically forbidden region in the (my 2 , mo)-plane. (With- 
out coannihilation, the requirement of the relic abundance of the neutralino to be consistent 
with that allowed by the age of the Universe (Q x h 2 < 0(1)) often provides a stringent up- 
per bound on the parameters m^ 2 and uiq (9|, [ll]].) The effect has since been included 
in a number of recent analyses, e.g. in [15, 16], and in a publically available package for 
computing the relic density micrOMEGAs [17]. 

One should mention that there are two other ways of evading this otherwise generic 
cosmological bound on and wiq. One is realized for mo 3> m^- In this region 

one invariably finds it difficult to satisfy the conditions of radiative electroweak symmetry 
breaking (EWSB); in other words the square of the Higgs/higgsino mass parameter \x comes 
out to be negative. In a very narrow corridor along the region of no-EWSB, /j, grows rapidly 
from zero but it is still less than m 1 / 2 [H, 18]. As a result, the LSP has a sizeable higgsino 
component (although it is still mostly a bino, like in the rest of the (mi/2, mo)-plane) and 
its relic density is typically small. In fact, because of the growing LSP mass and its gaugino 
fraction, the relic density increases rapidly along a very steep slope from very small values, 
characteristic for lighther neutralinos with a larger higgsino admixture, to larger (and often 
cosmologically excluded) values characteristic of heavier and bino-dominated neutralino. 
As a result, the cosmologically favored range 0.1 < £l x h 2 < 0.2 is only realized there for a 
rather narrow range of m-i/2 HI) 19]- 

The other important escape route from the cosmological bound on m^ 2 and mo occurs 
when tan (3, which is the usual ratio of the vacuum expectation values of the neutral Higgs 
scalars, is large, > 50. This is because the physical masses of the heavy Higgs scalar H and 
pseudoscalar A decrease with increasing tan (5 || . When the neutralino mass becomes large 
enough, close to half of the heavy Higgs boson mass, the LSP relic abundance becomes 
efficiently reduced through a relatively wide resonance involving mostly the pseudoscalar 
exchange. The effect is amplified by the coupling of A to down-type fermions which 
grows like tan/3. (The heavy scalar Higgs coupling also grows in a similar fashion but its 
dominant contribution is only p-wave and therefore suppressed by square of the WIMP 
relative velocity.) In the framework of the CMSSM, the effect is to open up wide fractions 
of the otherwise cosmologically excluded ranges of the (m^2, mo)~plane along the wide 
^-resonance (IJ, m, |(| ||, ||. The precise position of the resonance shows a sizeable 
dependence on some input parameters, most notably on the ratio mt/rrib, Ao, etc, and 
remains a subject of some ongoing debate. The full two-loop Higgs effective potential 
would have to be computed and implemented in the analysis to reduce the sensitivity to, 
e.g., the scale dependence. 

In the CMSSM a set of reasonably well-motivated unification assumptions leads to 
only four parameters: a common gaugino mass va\j2, a common scalar mass mo, a trilinear 
coupling Aq, as well as tan (3. One is also free to choose the sign of the //-parameter, while 
its magnitude is determined by the mechanism of electroweak radiative symmetry breaking 
(EWSB). The CMSSM can therefore be considered as a well-motivated SUSY model with 
the smallest number of independent parameters. (In particular, the so-called minimal 
supergravity (mSUGRA) model can be viewed as a specific realization of the framework.) 
The CMSSM has become a benchmark model for the LHC and other SUSY searches. 
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As mentioned above, the neutralino-stau coannihilation effect is particularly important 
in the framework of the CMSSM, but it can affect the neutralino relic density also in the 
more general MSSM. This will be the framework in which we will work here for the sake 
of generality. 

In this paper, we will present a full set of exact, analytic expressions for the tree- 
level cross section of the neutralino coannihilation with sleptons into all two-body final 
states in the general MSSM. In our analysis we will make no simplifying assumptions 
about the neutralino, nor will we assume the degeneracy of the left- and right-sfermion 
masses. We will not consider here the possibility of CP and flavor violation in the slepton 
sector although we will assume a general form of the left-right slepton mixing within each 
generation. We will include all tree-level final states and all intermediate states. We will 
also keep finite widths in s-channel resonances. A set of expressions for the neutralino- 
slepton coannihilaton was given in [14, 23 1 but only in the approximation of the partial 
wave expansion. Furthermore, these formulae did not include the effects of the tau Yukawa, 
of the T\ — T2 mixing and in some channels of the mass of the r, etc, which make them less 
reliable at large tan/3 > 20 |24]. 

The results presented here are exact, include all the terms and are valid both small 
and large values of tan (3, and both near and away from resonances and thresholds for new 
final states. This paper is meant to be a follow-up to ji where we have calculated all the 
analytic cross sections of all tree-level processes for the neutralino pair-annihilation into 
all two-body final states. We follow the same conventions and notations as in [jj]]. 

The plan of the paper is as follows. In Sect. ^ we briefly review the formalism for 
computing the relic density in the presence of both annihilation and coannihilation. In 
Sect. p| we introduce the relevant ingredients of the MSSM and list all the neutralino pair- 
annihilation channels. Explicit expressions for the coannihilation cross secion are given in 
Sect. |4|. In Sect. |5] we present some numerical examples and in Sect. || we summarize our 
work. Appendix A contains a complete list of relevant couplings and in Appendix B we 
provide expressions for several auxiliary functions used in the text. 



2. Calculation of the Relic Density with Coannihilation 

As summarized, e.g., in jjj, the time evolution and subsequent freeze-out of a stable relic 
in an expanding Universe are described by the Boltzmann equation 

*L = -3Hn - (av M0l ) [n 2 - (n^) 2 ] , (2.1) 

where n stands for the species' number density, n eq is the number density that it would 
have had if it had remained in thermal equilibrium, H{T) is the Hubble expansion rate, 
a denotes the cross section of the species pair-annihilation into all allowed final states, 
i>M0i is the so-called M0ller velocity |25| which is effectively the relative velocity of the two 
initial-state non-relativistic particles in the CM frame, and (<r«M0i) represents the thermal 
average of ctvmoI- hi the early Universe, the species were initially in thermal equilibrium, 
n = n eq . The number density decreased with the expanding Universe and the relic froze out 
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of the thermal equilibrium when its typical interaction rate became less than the Hubble 
parameter. 

In a more general framework with coannihilation, the stable particle x is the lightest 
of N species Xii each with mass mi, number density rij, equilibrium number density 
and the number of internal degrees of freedom g^. Griest and Seckel || showed that the 
total number density of all the species taking part in the coannihilation process 

N 

n = ^2m (2.2) 

i 

obeys the Boltzmann equation as given in eq. The quantity (at>M0i) now stands for 

eq n eq 

(<™ M0 i> = ^(ff.i%>^^, (2-3) 
ij 

where 



c-ij = <?{XiXj -> all) 



(2.4) 



and Vij are the relative (M0ller) velocities of the coannihilating particles. Edsjo and Gon- 
dolo recast eq. (|2.3|) into a convenient form involving only one-dimensional integrals 



(<™ M0 i) 



A, 



(2.5) 



where s = (pi + Pj) 2 is the usual Mandelstam variable, -fQ denotes the modified Bessel 
function of order i, and the kinematic factor (3f(s,mi,mj) is given by 



(3f(s,mi,mj) 



(mi +rrij) 2 



1/2 



(rrii — mjY 



21 1/2 



(2.6) 



Note that eq. ( p.6[ ) reduces to a familiar expression derived by Gondolo and Gelmini [25] 
in the case of single particle annihilation. (Compare, e.g., eq. (2.3) in Q.) Note also that, 
because of the assumed -R-parity, after the freeze-out only the LSP x wm survive and its 
number density will at the end be n x = n. 

Eqs. (^) and ( |2.5D can be rather accurately solved by iteratively solving the equation 

m x g cS 



In 



2tt 3 2 V 2g*G N 



45 1/2 
■{<rvu0i){xf)Xf 



(2.7) 



where g e R is the effective number of degrees of freedom of the coannihilating particles ||, 
g* represents the effective number of degrees of freedom at freeze-out {-y/gi ~ 9), Gn is the 
gravitational constant, x = T/m x and the freeze-out point Xf = Tf/m x is roughly 1/25 to 
1/20. (See, e.g., ' 
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The neutralino LSP relic abundance is £l x h 2 



p x j Pcriu where the critical density is 



Pcrit 



1.05 x 10 (h ) GeV/cm , and p x today is to a good approximation given by 

3 



Px 



1.66 



Mpi \T, 



^ ) T 3 ,/c, 



dx(av M 0\)(x) 



(2- 



where Mp\ = denotes the Planck mass, T x and T 7 are the present temperatures 

of the neutralino and the photon, respectively. The suppression factor (T x /T 7 ) 3 1/20 
follows from entropy conservation in a comoving volume [27]. 



3. WIMP Coannihilation in the MSSM 

We will be working in the framework of the general MSSM. (For a review, see, e.g., p8| . 
We follow the conventions of Iffi] .) We will not define it here but instead refer the reader 
to our previous publication Q where also all the relevant quantities are introduced. All 
the remaining couplings that we will need here are summarized in Appendix A. 

In this Section we summarize all the coannihilation channels of the neutralino LSP 
with the sleptons into tree-level two-body final states. They are listed in Table |l[ We only 
neglect slepton coannihilation with the lightest chargino and next-to-lightest neutralino 
which might be of some importance in the higgsino-dominant case. However, in this region 
the neutralino-chargino coannihilation is very effective anyway in reducing the relic density 
below any interesting level. 

We further neglect (s-)lepton flavor (generation) mixing, but include the left-right 
mixing for the sleptons. Thus a = 1,2, where the index a denotes the slepton mass state 
within each family, t and £' represent charged leptons of different generations. 



4. Exact Expressions 

We now move on to present a full set of exact, analytic expressions for the total cross 
sections for the neutralino-slepton coannihilation processes. We have included all con- 
tributing diagrams as well as all interference terms and kept finite widths of all s-channel 
resonances. In addition to neglecting CP violating phases in SUSY parameters we have 
also assumed no mixing among different generations of leptons and sleptons, although we 
have kept the left -right slepton mixing within each generation. The results presented here 
are meant to supplement and extend the ones presented in our previous paper Q for the 
case of the neutralino-pair annihilation. 

In the presence of coannihilation the formalism used in Q must be generalized to 
include different initial states. First, following (?], 23, [?(]] we introduce a Lorentz-invariant 
function Wij(s) 



Wij ( s ) = lf dLIPS \A(ij -» all)| 2 (4.1) 



where \A(ij — > all) | 2 denotes the absolute square of the reduced matrix element for the 
annihilation of two initial particles ij into all allowed final states, averaged over initial 
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spins and summed over final spins. The function vuij(s) is related to the coannihilation 
cross section aij(s) in eq. ( fEB] ) via |7|, [U|] 

w ij( s ) = -f3f{s,m u mj)aij(s). (4.2) 

The expression fl4.2|) readily reduces to a more familiar form in the case of single-particle 
annihilation. (Compare, e.g., eq. (4.2) of §.) 

Since Wij(s) receives contributions from all the kinematically allowed annihilation pro- 
cess ij — » /1/2, it can be written as 



W M = JhYl ce i s -( m h + m / 2 ) 2 ) Pf( s i m fi, m h)™ij->fih( s ) 



32 7T 

hh 



(4.3) 



where the summation extends over all possible two-body final states /1/2, m f 1 and mj 2 
denote their respective masses, and 

c = f c f if / 1(2) = /(/") 
1 1 otherwise, 

where cj is the color factor of SM fermions (cj = 3 for quarks and Cf = 1 for leptons). 

Since u>jj are Lorentz-invariant functions, we choose, for convenience, the CM frame 
in which the function w^^^f^s) can be expressed as 

1 f +1 

*%^/i/ 2 ( s ) = 2 J dcos9 C M \A(ij -> /1/2) | 2 , (4.5) 

where 9cm denotes the scattering angle in the CM frame. In other words, we write \A(ij — > 
/i/2)| 2 as a function of s and cos 9cm, which greatly simplifies the computation. 

We will follow Table |l| in presenting explicit expressions for Wij^f^f 2 (s) for all the 
two-body final states. All the couplings are defined in Appendix A. All other auxiliary 
functions, are listed in Appendix B. Some symbols are obvious {e.g., my/ is the mass of 
the W^-boson) and will not be defined here. 



We begin by presenting the results for l a l* h — > vector boson pairs: WW, ZZ, Zj and 



77. 

4.1 l a l% 



WW 



This process involves the s-channel CP-even Higgs boson (h and H) exchange, the four- 
point interaction, the s-channel Z-boson and photon exchange, and the u-channel sneu- 
trino {yi) exchange 



~__ _ ~(h,H,P) ~{Z,i) ~{u) ~(h,H,P-V) ~(Zd-v) 

w i a q^ WW — W WW ~r W WW ~r W WW W WW " "' > 



J WW 



(4.6) 



Higgs (h,H) exchange (+ Point interaction): 



w 



(h,H,P) 
WW 



r tWWrr , tU a r 

(jtbtuM, w 

r=h,H 



s — ml + iT r m r 



5 2 — \mLrS + 12m 4 



w 



4m 4 , 



(4.7) 
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• Z, 7 exchange: 



71' 



(2,7) 



S — 777-7 + 



+ 



{s-(m~ +m~) } 



. 2 s 3 + lQm 2 w s 2 - 68m^s - 48tt7^ 
x{s - (mj - m~ ) } • n 



(4.8) 



sneutrino (yi) exchange: 



w 



(V) 
WW 



777 



w 



c qv e w- c v*e a w+ 



X [T4 - 2(m| + 777| + 2777^)T 3 

+{7771 + m| + 4mf o m| + 2m 2 w (mf a + m| + 3t77^)}t 2 
-2{(m| + m|)(m| a m| - 777^) 
+m w (m~^ + 777- - Am-^m-^ + 2777^) j 71 
+(777^ -m 2 w ) 2 K -"^ol, 



(4.9) 



where the auxiliary functions 7$ (i = 0, 1, 2, 3, 4) are defined in Appendix B. The arguments 
for these functions in the above should be understood as % = %(s, m~ , m- , m^, m 2 ^, 
m~ , m~). Throughout the text, the second and the third arguments denote the masses 
for the two initial particles. The fourth and the fifth ones correspond to the masses for 
the two final particles. The last two arguments are the masses for the exchanged particles 
in the relevant t- and/or n-channel diagrams. Assuming this convention, we will omit the 
arguments for % in the following. 

• Higgs (h, H) (+ Point) - sneutrino (pt) interference: 



w 



(h,H,P-v) 
WW 



2777^ 



Rc 



E 



(jWWr(ji* b i a r 

s — m 2 + iY r m r 



c e* b e a ww c tiv t w- C v*iaW^ 



2 1 2 

2 2 \ / 2 2 

777- — 777w ( 777- — 777- — 777 



S + S(m~ + 777- - 2777-^ - 4m w ) + 2777^/(777- + 777- - 2777-^ + 2777 

-2[s(m 2 7 



w) 



AMI 2. 

m w)( m j b - m l 



w) 



+2m 2 v {m^ + - m\ mj - {m\ f + m 2 w )(m\ + m 2 ~) + (m 2 ~ i - m 2 w ) 2 }\T 



(4.10) 



where the auxiliary function T is defined in Appendix B. The arguments for these functions 
in the above should be understood as T = F(s, m~ , m~ , m"^, m 2 ^, m~ ). Throughout the 
text, the second and the third arguments of T denote the masses for the two initial particles. 
The fourth and the fifth ones correspond to the masses for the two final particles. The 
last argument is the mass for the exchanged particles in the relevant t- and/or Ti-channel 
diagrams. With this understanding, we will omit the arguments for T in the following. 
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• Z, 7 - sneutrino (i/g) interference: 



w 



WW 



6m^ 



Re 



(jWWZ(jtlt a Z 
s — m? z + iTzmz 



+ 



e 2 5 ab 



c e* b u e w~ c v*e a w+ 



x[l^ w + 6I^ w T], 



(4.11) 



where 



tWW 



tWW 
1 2 



-s 3 + 2s 2 (m~ +m~ 
+s{-m~ 



9m^) 



4 

m 4 



lOmimi + 12(m£ + 2m^)(m£ + m| 



+4(-3m~ + §m w m~ i + 7m^/)} 
+2m 2 w {-9(mi + mi ) + 30m 2 ; mf 

-4(3m§, + m^)(mf + mf ■ ) + 12(m|, - m^) 2 } 



-- m w( m j 

s la 



2 \2 

m 4 } ' 



2/2 2 2\/2 2 2 \ 

-s (m~ - m~ - m w )(m~ b - m~ - m w ) 



,2 2 



+s{m| j m|(m| a + m|) - (m| + 3m w )(m^ + m| ) 
-2(2m| £ + m w )m\m\^ + (3mi + 5m^)(m£ + m|) 
-2(^| - 4m w mj, e + m^m 2 ^ + 2m^)} 
+2m 2 y{(m| i + m~ )(m~ + m~ - 3m~m~) 
-(m§ £ + m w )(mj a + m|) + 4m|m|(2m^ + m 2 ^) 
-(m| + m| )(m§, - m 2 y )(3m| / + m 2 ^) + 2(m^ - m 2 ^) 3 }. 



(4.12) 



(4.13) 



4.2 £ a £* 



This process proceeds via the s-channel CP-even Higgs boson (h and H) exchange, the 
four-point interaction, and the t- and it-channel slepton (£ a , a = 1,2) exchange 



*zz ~ w zz 



ir, .,- - w-,; ' + w zz 
• Higgs (h,H) exchange (+ Point interaction): 



(t) ~(h,H,p-e) 



(4.14) 



71' 



(h,/r,p) 
zz 



C ZZr C £ b e " r 



r=h,H 



s — m 2 . + iT r m r 



tUaZZ 



s — \mr z s + 12m] 
8m| 



(4.15) 



slepton (£ c ) exchange: 



w 



ZZ 



mi 4^ 



c,d=l 
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T 4 - 2(m| + ml + 2m|)T 3 



+ {mi + mi +4mjmf + 2m z (mj + mi + 3m|)]T 2 

-2[(mf +mf)(mf mf -m|)+m|(m| + mf - 4m? m? + 2m|)]7i 



+ (mi -m|) 2 (m| - m|) 2 T 

,2 i _2 



34 + [s(mj + m£ - 2m|) - 2(m£ - m|)(m| - m|)]^ 2 



s 2 (m~ -m z )(mf - m z z ) 



+s{-m~ m~ [m~ +m~) + 3>m z {m~ + m~ ) - 3m z (m 7 + m~ ) + 2m z } 



+« -m|) 2 (mf -m|)^ 



2 n2i 



(4.16) 



Higgs (h, H) (+ Point) - slepton (£ c ) interference: 



w 



(h,H,p-e) 



zz 



2m% - 

z c=l 



(-iZZrcilUar 

E ° ° 6 Ql*l a ZZ 

s — ' n 



,r=h,H 



s — ml + iT r m r 



c e* b e c z c e*j a z 



x \ s 2 + s(mi +mi -2mi -4m z ) + 2m z (mj + mi -2mi + 2m|) 
-2[s(mj -mj -m 2 z )(mj - mj - m|) 
+2m 2 z {mj + m~ — ra~ ra~ + ra~ (m~ — m~ — m~ — 2m^) 



-m\{mi + m 2 ~ - m\)}\T 



(4.17) 



4.3 ^a^-fr 



Z7 



This process proceeds via the four— point interaction, and the t— and u— channel slepton (t a 
a = 1, 2) exchange 



~(P) ~(^) ~(P-<?) 



(4.18) 



• Point interaction: 



4? = 3 



(4.19) 



slepton (£ c ) exchange: 



w w = JL v c r ° taZ c e * e * l c e y aZ *c e i t ' n * 

^7 m 2_ 



z c,d=l 



> j -T 3 * + (2m| - mj b + 2m|)^ 



_ {(m 2 _ m | ) 2_ 2 2 (m 2 +m | )} ^_ m 2 (m 2 _ m | 



2 



2 



Z* 



m . 

^ c,d=l 
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-T 3 U + (2m| - m| + 2m|)T 2 " 
-{(m| - m|) 2 - 2mf a (m| + m|)}7J« - mf a (m| - m|) 2 T « 

+J_ Re y c i ^ z C I ^C I y^*C I ^ z *( S -2m 2 7 - 2mj + m|) 

z c,d=l 

-23>2 + -(a + m|)(ro2 -mj)y 1 

+ [5(771- + 777- - 2m|) - (777- + 777- ) + 3m|(m~ + 777- ) 



(4.20) 



Point - slepton (£ c ) interference: 



rXP-t) _ 1 



■77? 



Z c=l 



3t77j + 777j + 2m~ — 2m z m z (m~ — m~ ) 

*-a 4 4 S 7a 4 



+ [s(m~ -m~ + m 2 z ) + (2m~ - m z )(m~ - m 2 z ) 



+mf (2ml - 4m| - 3m|) - 2m 2 z mj\P 



+ _L Re ^4.^*^4.7^42 
m 2 z ^ 

3777- + 777- + 2777- — 2m| H m|(m~ — 77l~ ) 

4 *-a 7 C 5 "a ' ' 

2 ™2 , 2 \ , 7o™2 2\/ 2 2 



X 



1 



+ [s(m| - m| + m|) + (2m| - m|)(m| - m|) 

+m| (2m! - 4m| - 3mi) - 2m 2 ,m| 1^* 

4 V 4 4 z/ z 4 J 



(4.21) 



4.4 £ a £* b — 77 

This process proceeds via the four— point interaction, and the t— and it— channel slepton (i-ai 
a = 1, 2) exchange 

= «;£> + eg + : 



77'; 



+77 77 1 77 ' 77 



Point interaction: 



slepton (t a ) exchange: 



w 



(P) 

77 



8e 5 ab ; 



• Point - slepton (£ c ) interference: 



4(T 2 + 2mf Ti + mf T ) - (s - 4m^ )0>o 



^ 7 - £ ) = 2e 4 5, 



ab 



-4 + (s - 8m~ )J C " 



(4.22) 
(4.23) 

(4.24) 
(4.25) 
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Next, there are processes of the type l a l\ — > vector boson - Higgs boson: W ± H T , Zh, 
ZH, ZA, 7/1, jH and jA. 

4.5 lj* h -> W±ifT 

The process 4^ — ► W + H~ involves the s-channel CP-even Higgs boson (h and H) and 
CP-odd Higgs boson (A) exchange, and the u-channel sneutrino (ug) exchange 



*W+H- 



~~ "V+h- ^ "V+h- ^ "V+h- 



(4.26) 



Higgs (h, H,A) exchange: 



w 



(h,H,A) 
W+H- 



c w-H+r C e* b e a r 

2^ o _ ^,2 



s — m 2 + iF r m r 



s 2 - 2s(m 2 H± + m 2 w ) + (m 2 H± - m 2 w ) 2 



m 



W 



(4.27) 



sneutrino (vg) exchange: 



~{y) 
W W+H~ 



rn 



w 



c v*e a w+ c e* b p e H- 



T - - 2(m| + m%y)7? + (m| - m 2 w ) 2 T « 



(4.28) 



• Higgs (/1, i?, A) - sneutrino (ug) interference: 

(jW-H+i 



~{h,H,A-v) fL Rc 

W W+H~ ~ ™2 



rn 



w 



,r=h,H,A 



c qe a r 



s — m 2 + iT r m r 



2 1 2 . r / 2 , 2 2 

s - m H ± + m^y + |s(m^ + m w - m- 



+m 2 w {m\ - 2m| + m 2 ~ i - m 2 w ) 



+m 2 H± (m 2 ~-m 2 ~ i + m 2 w )}T u 



(4.29) 



The contribution from £ a £^ — ► H^ can be obtained by interchanging the indices a and 
b, but this is obviously equal to the one above: 



w e a e b -->w-H+ w i b i* a -*w+H-- 



(4.30) 



4.6 i a i* h — ► Z/i, Zii 

The process — > Z/i proceeds via the s-channel CP-odd Higgs boson (A) exchange, the 
s-channel Z-boson exchange, and the i- and ^-channel slepton (£ a , a = 1,2) exchange 

w~ ~ - w (A) + w (Z) + w {1) + w {A ~ Z) + w {A - J) + w {Z ~~ t] • (A 31) 

w g a g*^Zh — W Zh + w Zh + W Zh + w Zh + w Zh + • K^ 1 ) 

• A exchange: 



(A) 
Zh 



C 



ZhA 



QlllaA 



s — m A + iT avtia 



s 2 — 2s{w? h + m 2 z ) + {w? h — m 2 z ) 

171% 



(4.32) 
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Z exchange: 



w 



I2m% 



s — m z + iTzmz 



{3(m| a -m|) 2 +m|} 

-2s|3(m~ - m~ ) 2 (m| + 2m|) + m z (m 2 + m~ + m 2 - 5m|)| 

+m z {(m 2 h - m|) 2 +4(m| + m|)(m^ - 5m|)} 
+(mf -mf ) 2 (3mt + 6m|m^ + 19m|) 

|(m£ - m|) 2 (3m^ - 2m|m^ + m|) + m|(mj[ - m z f(m^ + m|)} 



2 2 

™Z 



S 

4 



+ 4^1(^-1 -m|) 2 (mf -m~) 



2 \2 



slepton (£ c ) exchange: 



(4.33) 



w 



mi 4-^ 



c,d=l 



x [7? - 2(mf o + m|)T/ + - m|) 2 r o * 
! 2 

_| t_ C e * £ah C £ b C l '<l lah *C l b id Z * 



z c,d=l 

x [r 2 « - 2(m| + m|)T« + (m| - m|) 2 T c 
! 2 

+ _R e V" c i*jaZ C ete c h c r d £ a h* c qe d z* 



x \-2y 2 + -(s-mi + m z ){m\ a - mjjy, 

{/I \ f A A \ 

s(rrij + m~ — 2m z ) — (m~ + m~ ) — 2m z m h 
ta *-a *-b 



+(3m| - m|)(m| + m|) + -(m^ - m|)(m| o - m|) 2 



•^(m 2 -m|) 2 (mf a -m|) 2 }% 



A - Z interference: 



(A-Z) 
Zh 



srrir 



-Re 



C-ZM^/aA \ * / c zzh C r b taZ \ 
s - m 2 A + iTAiriA J \s — m z + iF z mzJ 



x(mj - mj )(s - m 2 z ) [s 2 - 2s(m 2 h + m|) + (m 2 h - m|) 2 ] ; 



(4.34) 



(4.35) 
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Higgs (A) - slepton (£ c ) interference: 



w 



{A-t) 
Zh 



m, 



(jZhA(je;e a A \ ' 



^ ys -m 2 A + iT A m A J 



c r c I a z c ItI ch |_ (s _ m 2 + m 2 z) + [g(m | _ m | _ m 2 z) 

-m\{m\ + m|) + (m| - m\){m\ - m 2 z ) + 2m 2 z m^J rt ^ 
_ c e*J a h c qI c z _ m 2 + m | ^ + [ s ( m | -m~ - m|) 

-m|(m£ + m|) + (m| - m|)(m£ - m|) + 2m|m|j^ u } 
Zl - slepton (^ c ) interference: 



w 



(Z-t) 
Zh 



m. 



-Re£ 



2 / (jZZhQtlt a Z \ ' 



c=l 



s — m c 



c-w c ^;_ (s _ m 2 )(m 2 



+2m|(m~ - m~ + m|) - -m\{m\ - m 2 z )(m~ - 



+ 



s{m 7 —mj +m z )(m 7 — m 7 -m z ) 



m 7 



+m|{-2mj + rrij (3m~ + m~ + 3m 2 z ) — Srrij 

-f-c *-a -t-Jj ""' 

+3m~ m~ — 2m~ + m^(4m~ + m| ) — m^} 
-m 2 h (mj a - m| + m|)(m| - m| + m|)] J* 7 *} 
+ ^a# z f( S -m^)(m? -m?) 



1 



2 \ i 2/2 2 

z) + - m zi m h ~ m z 



+2m z (m~ -m~ +m z )-\ — m z (m h — m z )(m~ — m~ 



+ 



+m z {-2m~ + m~ (3m~ + m~ + 3m|) - 3m~ 
+3m~ m~ - 2m~ + m|(4m~ + m~ ) — m z } 

+m 2 h (mj a - m| - m|)(m| - m| + m|)] J^} 



The expression for w^r j*_ >ZH can be obtained by a simple substitution h — ► if. 



4.7 t a (.^ 



ZA 



(4.36) 



(4.37) 



This process proceeds via the s-channel CP-even Higgs boson (/i and if) exchange, and 
the t- and n-channel slepton (£ a , a = 1, 2) exchange 



(h,ff) 



(h,H-l) 



(4.38) 
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• Higgs (h, H) exchange: 



w 



(h,H) 
ZA 



C ZAr C Ife a r 



. rT s — mi + iT r m T 

r=h,H 



s 2 — 2s(m 2 A + m' z ) + (m A — m 2 z ) 2 

2 

m% 



(4.39) 



The expression for w ZA can be obtained from w Zh by a simple substitution h ^ A. 
• Higgs (/i, H) - slepton (l c ) interference: 



w 



(h,H-t) _ 2 



ZA 



^ c=l \r=h,H ' 



C lfi a Z c Ifi c A (_ {s _ m 2 A + m 2 z) + [s{m 2 _ m 2 ^ 

-m~ (m^ + m|) + (m~ - m 2 z )(m A - m|) + 2m|m|]J r *| 
^e&A^z f_ (s _ m 2 + m | } + [s(m 2 _ 2 m | } 

-m\{m\ + m|) + (m| - m%){m 2 A - m|) + Im^m 2 -]^}^ . (4.40) 



4.8 ^ a €* -> 7/1,7// 

The process — »■ 7/1 proceeds only via the i- and u-channel slepton (£ a , a = 1,2) 
exchange 



(4.41) 



slepton (£ c ) exchange: 



2 o™2 



(7? + mf^) + (7? + mf T ") 



+(s + m£-2mf -2m^)3; 



(4.42) 



The expression for w~ c tt* „ can be obtained by a simple substitution h — > H. 
4.9 IJl -> 7 A 

This process proceeds only via the t- and u-channel slepton (£ a , a = 1, 2) exchange 

_ ~(e) 

w^ A can be obtained from wfy by a simple substitution h — > A. 



(4.43) 



Next we proceed to present the results for i a l* h — > Higgs-Higgs pairs: 
/tA tfA, AA and iJ+i/-. 
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4.10 l a l% -» hh, HH, hH 

The process £ a £t — > proceeds via the s-channel CP-even Higgs boson (/i and i?) 
exchange, the four-point interaction, and the i- and u-channel slepton (£ a , a = 1,2) 
exchange 



~(h,H,P) . . ~(h,H,P-t) 

W e a l* b ->hH = W hH + W hH + W hH ■ 



Higgs (h,H) exchange (+ Point interaction): 



w 



(h,H,P) 



hH 



(jhHr (jt%tar 



rr s — mt: + iT r m r 

r=h,H 



eUahH 



slepton (£ c ) exchange: 



(4.44) 



(4.45) 



w 



(!) 

hH 



c etiah c eie c H c e* d e a h* c e* b e d H*j-t + ^ c e* c e a H c efe c h c e d e a H* c r b e. d h*j-u 



c,d=l 



c,d=l 



-2Re C e * eah C r b e cH C r d e a H* c e* b e d h*y o . 

c,d=l 

Higgs (h, H) (+ Point) - slepton (£ c ) interference: 



(4.46) 



c=l \r=h,H 



(jhHr (j£ b £ a f 

s — ml + iT r m r 



UahH 



^(je.* c iah(ji* b icHj:t _j_ QV c t a H (ji* b tch -pu 



(4.47) 



The expressions for /i/i final state are obtained from the above by replacing C hHr , 
C r b LH^ c e*J a H^ C hH~t~t a with C rfthr ) c,?;?c/» ) C hft? 6 ? », respectively, and multiplying 

w by a factor of 1/2 for identical particles in the final state. The contributions for HH 
final state are obtained in an analogous way. 

4.11 IJ* -» hA, HA 

The process l a Pt — > /iA proceeds via the s-channel CP-odd Higgs boson (.A) exchange, the 
s-channel Z-boson exchange, and the t- and ^-channel slepton (£ a , a = 1,2) exchange 



^7 7. 



<„«-M ~ W M + ^ + w ^ + w hA 



+ w hA + w 



hA 



(4.48) 



A exchange: 



w 



(A) 
hA 



(jhAAQl* b laA 



s — m A + iTAfnA 



(4.49) 
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Z exchange: 



w 



(Z) = 1 
hA 3s 2 m% 



(jhAZ Qt* h l a Z 



s — m 2 z + iTzmz 
x \s 4 m z — 2s i m z (m~ + m~ +m\ + mf l ) 
+s 2 {(ml - m|) 2 [3(mi - m 2 ) 2 + m%\ 
+Am A z (m 2 A + m 2 h )(mj a + m|) + m|(m^ - m 2 ,) 2 } 
-2sm| |(m| - m|) 2 [3(m^ - m£) 2 + m|(m^ + m 2 h )] 
+m 2 z (m 2 A - m\f(m\ a + m|)} 
+4ml(m 2 4 - m 2 ) 2 (m| - m|) 2 ] ; 



(4.50) 



• slepton (^ c ) exchange: 



~W = c e* c e a h c r b e c A c e* d e a h* c eie d A*^t , >p c r c iaA c e* b e c h c e*e a A* c e* b e d h*q-u 

hA / j ' / ^ 



c,d=l 



c,d=l 



-2Re ^ c^^c^c^ A *cW*;); ; 

c,d=l 



(4.51) 



A - Z interference: 



w 



{A-Z) __ 
hA 



2Re 



(jhAA C e*e a A \* I C hAz c qe a z \ 
s - m\ + iTatha J Is- m 2 z + ir^m^ i 



(m 2 A - m 2 h )(mj - mj )(s - m|) 



(4.52) 



• Higgs (^4) - slepton (£ c ) interference: 



w 



(A-e) 

hA 



C hAA C r b e a A \ : 



^ \ s - m\ + iTAmA I 
Z - slepton (£ c ) interference: 



(4.53) 



(Z- 
hA 



c=l 



QhAZQl* b l a Z 



-V 



(jt* c tah(jt b tcA 



s — m 2 z + iTzmz I 
x [2m| - {(mi - ml)(mj a - m|) + m|(mf a + m|) 
+m 2 z (-s - 2m 2 ~ +m 2 A + m 2 h )}F^ 
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mi ^— ' 



Z c=l 



(jhAZ Qt^l a Z 

iT z m z J 



s — mi, + il 



x ^2m| - {-(m^ - m 2 h )(m^ - m|) + m|(n^ + m|) 
+m|(-a - 2m^ + + m|)}^ u . 



(4.54) 



The expression for wj j*_^ HA can be obtained by a simple substitution h — > iT. 



4.12 £ a £* -> AA 



This process proceeds via the s-channel CP-even Higgs boson (h and -fT) exchange, the 
four-point interaction, and the t- and u-channel slepton (£ a , a = 1,2) exchange 



~{h,H,P) ~{t) ~{h,H,P-t) 
w Ul -^AA = W AA + W AA + W AA 



(4.55) 



Higgs (h,H) exchange (+ Point interaction): 



w 



(h,H,P) _ 
AA - 



r<AAr riRlar 

^2 —4 c e * £aA - 



r=h,H 



s — mi + iV r m r 



(4.56) 



slepton (£ c ) exchange: 



w 



(4.57) 



c,d=l 



Higgs (/i, i?) (+ Point) - slepton (£ c ) interference: 



c=l 



>r=h,H 



s — + iT r m r 



c £*£ a A c eie c A 



T. 



(4.58) 



4.13 i a i* h -> Jf+J/ 



This process proceeds via the s-channel CP-even Higgs boson (h and H) exchange, the 
four-point interaction, the s-channel Z-boson and photon exchange, and the u-channel 
sneutrino {yi) exchange 



w 



Ah,H,P) 



-XZ,y) 



w H+H J + + w)j' +H _ + w H+H 



{h,H,P-u) + ~{Z„-v) . 



H+H- 



(4.59) 



Higgs (h,H) exchange (+ Point interaction): 



w 



(h,H,P) 
H+H- 



_^ ° b _ c £* b £ a H+H- 

r=h,H 



s — mi. + iT r m r 



(4.60) 
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Z, 7 exchange: 



w 



H+H- 



1 

3s 



C H+H-z c e* b i a z 2^ 



s — m z + iFzmz s 



x(s - 4m 2 H ±)[s - (to- + m~J 2 ][s - (to- - to- 



sneutrino exchange: 



w 



H+H- 



(jv*t a H+ c r b u e H- 



Higgs (h, H) (+ Point) - sneutrino (ug) interference: 



(4.61) 



(4.62) 



wf + H H P ~~ v) = 2Re 



r H+H-r r t*£ a r 

q . " 

,r=h,H ' 



s — m 2 . + iT r m r 



(4.63) 



7 - sneutrino (z^) interference: 



= 2Re 



C H+H-z c r b e a z e 2 S 



+ 



— m? 7 + ir^m^ s 



2 + (s + 2m§ £ - 2to^± - mf - mf ).F 



(4.64) 



Finally, we present the results for i a l\ — > // where / denotes any of the SM fermions. 
4.14 lj* — > 

This process involves the s-channel CP-even Higgs boson (/i and H) and CP-odd Higgs 
boson (A) exchange, the s-channel Z-boson and photon exchange, and the i-channel 
neutralino i = 1,2,3,4) exchange 



~(h,H,A) . ~(Z,i) . ~(x°) . . ~(h,H,A- x °) , ~(Z, 7 - X °) /, «^ 



Higgs (h, H, A) exchange: 



w 



(h,H,A) = 2 



7 exchange: 



C Ur C ttUr 



, s — m| + iT r m r 

r=h,H 



(s - Am}) + 2 



s - m\ + ir^m/i 



(4.66) 



w 



C uz c r b t a z eeeSab 



s — m\ + iTzmz s 
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x (s + 2m\ ) { s - (m~ + m~ ) 2 } J 1 - 



(m~ 



3sm z 



C UZ C tttaZ 



s — m| + iFzmz 
s 3 m z — 2s 2 {m z (m~ + m~ + 2mf) — 3mj(m~ —rn~) 2 } 



+sm z {{m z - 12m e )(m^ - mjj + 8m e m z (m~ a + m~ )} 



-\-2m i m z {m~ — m~ 

*~a *-b 



2 \2 



(4.67) 



where = — e < is the electric charge for the charged lepton t in the final state. For 
other fermion pair production, this should be replaced with the electric charge of the final 
particle. 

• neutralino (%°) exchange: 



w 



(x°) 



-2m\ [cl b R J LR m x om x o% + C^Wi] + Cf^m^m^s - 2mj)% 
'cSLlVi - (ml ~ mj)%] + C^ijT! - (m| - m, 2 )T } 



— m^m v o 
— m^m v o 

A 7 



abij 



(4.68) 



where 



^yabij 
LLLL 



C 



abij 
LRLR 

abij 
LLRR 



L L 

■xf^* r xf^ r xfV* r xftl 

L °R U L U R 



f~iabij 
LRRL 



c 



L L ^R ^R 

xftj^xftit^xf^^J 



+ (L^ R), 
R), 



°R U R 



c 



abij 
LRLL 



c xfrj* c ^Iit c x%t* c x%z + (L 



f-tabij 
RLLL 



c 



L ^R L L 

xftj* c xftlt c x%i* n x%t 



R 



Cr 



R), 

R), 
R), 



C 



abij 

LLRL ~ 



c 



x\t a i* nX \t b i n x%i* n x%t 



+ (L^ R), 



(-tabij 
LLLR 



C f^cf» e cf te *cf" e + (L~R); 



A — Z interference: 



w 



(A-Z) 



-8Re 



c UA C qe a A 



c uz c e h t a z 



s — m A + iY A"rfiA ) \s — m z + iT z mz 



m e (ml - m| fc )(a-m|); 
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(4.69) 



Higgs (h, H, A) - neutralino (%°) interference: 



w 



(h,H,A- x °) 



C Ur C lll a r 



i=l r =h,H 



s — raj. + iT r m r J 
C%m E [ - 2 + (ml + m| - 2m 2 „ - 2m 2 )^} + C%m x o(s - 4m 2 )F 

D%m e (mj -ml)+D%m^s 



4 / r UA r i1i a A \* 

+4Re^ _ p ~ I 



where 



Z, 7 - neutralino (x^ 1 ) interference: 



I s - + iTArriA ) 



c: 



(4.70) 

(4.71) 
(4.72) 



w 



c uz c r b I a z eeA 



ab 



4ReV , 

^ \ s - m z + iT z m z s J 

C%{ -s + ml+ m| - 2m 2 o - 2m, 2 

+2[-sm 2 c o - (m~ a - m 2 o)(m| - m 2 o) + m^j 

• 2m e m x o | - 2 + (s + m| + m| - 2m 2 + 2m^)^} 
^ / c uz c tii a z \* ! 

+ 4Re E ^2 ^ -2" 

\ s — m z + il / m z 

D%{m 2 z {-s + mj a + m| - 2mJ„ + 2m £ 2 ) 
-2[m 2 z {sm 2 x o + (m?r - m 2 ^)(m\ - m 2 )} 



7„ X^' 



+m 2 e {(mj -mj) 2 -m 2 z (mj + mj + 2m 2 -m 2 e )}}F\ 



-D% ■ 2m e m x o{mj a - m|)(s - m 2 z )T 



(4.73) 



4.15 l a i* h qq 

This process involves the s-channel CP-even Higgs boson (h and H) and CP-odd Higgs 
boson (A) exchange, and the s-channel Z-boson and photon exchange 

Ah,H,A) ~(z l7 ) , ~(A-Z) . 



W 



qq 



+ Wqq + Wqq 



(4.74) 
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Each contribution can be found from the respective ones for the li final state by replacing 
I with q. 



4.16 l a ll 



This process involves the s-channel Z-boson, and the t-channel chargino (xf, k = 1,2) 
exchange 



(4.75) 



Each contribution can be found from the respective ones for the it final state by replacing 
I with i/£. 



4.17 £ a £'* 



This process involves the t-channel neutralino (x®, i = 1,2,3,4) exchange 



1"2 



(4.76) 



neutralino (xf) exchange: 



w 



(x°) 
w 



where 



C 



abij 
LLLL 



- T 2 * -(s- m\ - mf )Tl - (mf - m 2 e )(rrA - m|)T * 



—2m£in£r 



C LRLR m x* m X^ + C TlRR T 1 + C LRRL m x'l m x ( ] ( S ~ m i ~ m D T 



abij 



-m v o 



- m x° 



Ct% L mATi ~ (mi - mj)Tt} + Cg* £ ro,{7? - (mf - m, 2 ,)^} 



abij 



(4.77) 



(~iabij 
LLLL 


(jXi Z a Z* (J^i ^ 


c x%z* c x%*t' 


+ 








T^abij 
LRLR 


r xft*j* r x^:z 

— °L U R 


c xf*J* c x%*t' 

L R 


+ 


(L 


<— > 




p^abij 
LLRR 




n xf*J* x %n> 

U R 


+ 


(L 


<— > 


fl), 


p^abij 
LRRL 




r xfV* r x%*t' 


+ 


(L 


<— > 


R), 


T^abij 
LRLL 


r x?rj* r x%*i 


c x°e*J* c x%*e' 


+ 


(L 


<— > 


R), 


T^abij 
RLLL 




-v°7*f* V /*/' 
(jXj*-a*- * (jXj*- h t 


+ 


(L 


<— > 


R), 


T^abij 
LLRL 






+ 


(L 


<— > 


R), 


T^abij 
LLLR 


= cf J * J *c^ :i 




+ 


(L 


<— > 


R). 
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4.18 i a t'* 



This process involves the ^-channel chargino (%^, k = 1, 2) exchange 



The expression for ) can be found from u>j;|, J by replacing £ and x° with vi and x^, 
respectively. 



-Xx°) 



(4.78) 



The second class of processes involves £ a £b — > and £ a £' b — > For each there is a corre- 
sponding process involving particles with the oposite electric charge. 



4» 19 ^'CL^'b ^ 

This 



process proceeds only via the i- and u-channel neutralino (x°, z = 1, 2, 3, 4) exchan 

w e a e b -*u ~ w u '■ 

• neutralino (%°) exchange: 



tge 
(4.79) 



~(x°) 
w~~ 
Lih- 



= E 



-2m? [^SW™*?™^ + D lrrl T i 



n abij 
' LRLR 



-T 2 -(s-mj -mj )Ti - (mj - m\)(mj - mj)T 



D a ^ LL {Ti - H a - m ^ + ^Su* - K - m?)To} 



+ 



- 2D L b LLL m x om X ( S ~ 2m e)y0 



2m £ 2 )3^o 



, 9 7->ab«j 
' LRRL 



-y 2 - {mj mj - mj)y 
-m e m x o Dfl 3 LR {y l + (s + m~ - m| - 4mf)>'o} 

+ (* " ™| + ™| " *rn 2 e )y Q } 
+m,m x o [l^i^i - (s + mf o - m| - 4mj)y } 

+D a R L j LL {yi ~(s- mf o + m| - 4m?)3; }^ 



(4.80) 
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where 



D 



abij 



D 



LLLL 

abij 
LLRR 



p.abij 
LRRL 



D 

LRLR 

n abij 
LLLR 



D 



abij 
LLRL 



j^abij 
LRLL 



D 



C 



^*cf^cf li cf J + {L^R), 



C 



L L 



^R U R 



+ (L^ R), 



C fV*cf^cf l 'cf* J + {L~R), 



abij _ c x^* r x ^* r x °i l b l r x °^ J 



C 



R 



c: 



abij _ 
RLLL — 



C£ a + (L R), 

C f^cf^cf^cf J + {L~Rl 
c fie* c fte* c ft* c f*J + {L ^ Rl 

C f^cf^cf V cf lt + {L~ R), 
c fte* c ftt* c ft* c f*J + {L ^ R) 



(4.81) 



4.20 £ a £' b -> it' 

This process proceeds only via the t-channel neutralino (x% i = 1,2,3,4) exchange 



(x°) . 



(4.82) 



neutralino (x®) exchange: 



w 



(x°) 



D 



abij i 2 2 \n~t 

LLLL m x° m x S S ~ m i- m t')% 



D LLRR m xl m X^ + D LRRL r i 



abij 



Xj X 



+D 



abij 
LRLR 



T 2 * -(s- mj - mj )T£ - (mj - m?)(m| - m £ 2 , )T * 



— m v o 



DfUnm^Tt - (mj - mj)T^} + D a ^ RL m e {^ - (mj - m%)7^ 



where 



— m v o 



Df^LL^in - (mj ~ mfK} + DZ J L L me{7? ~ (mj - mf,)^} 



\abij 



(4.83) 



fichu = c x°e* a e* c x°t b *e'* c x%e c 



D 
D 
D 



LLLL 

abij 
LLRR 

abij 
LRRL 

abij 
LRLR 



n abij 
LLLR 

r>abij 
LLRL 



r x°rj* r x°l*e'* r x /*j x °e b *t 

L/^ L D D 



' R 



'R 



r xXi* r x^e'* r x%£ r x%*i 



c 



c l a cr R 



r x\ilt* r x%*V* r x%i r x%* 

- ' 1 ^R 



c 



x\t a i*^ x %^'*„x%i^xfct 



+ (L^R), 
+ (L^ R), 
+ (L^R), 
+ (L^ R), 
+ (L^ R), 
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LRLL 
\abij 



d rlll = C^ PJ *C^ tb ' e *C X ^' J C X ^ Cl ' + (L <-► R). 



Finally, there are slepton-neutralino annihilation channels into a lepton and either a gauge 
or a Higgs boson. For each, there is an analogous one involving the opposite electric charge. 

4.21 I aX -> Z£ 

This process involves the s-channel lepton (£) exchange, the t-channel slepton (£ a , a = 1, 2) 
exchange, and the u-channel neutralino (%°, z = 1,2,3,4) exchange 



«^ oX _z/ = «>z< + ™zi + w z/ + W Z£ + w zt 



+ u> 



z^ 



lepton (I) exchange: 



w 



z ' " 2m|(s-m2)2 

x (5 — m~ + m?,){(s — raj) 2 + m^(s + m?) — 2m|} 



+4C^ a 11 Cf z m^m x {(s - mj) 2 + m z (s + mj) - 2m%] 
-12C™ u C e J z mjm 2 z (s - mj + m 2 ,) - 12Cl a u C e J z m e m x m 2 z (s + mj) 



where 



C: 
D 



uz 



uzf ± (c ^z )2j 

•ttZriUZ. 



2c\rc 



V U A » 



slepton (4) exchange: 



(4.84) 



,(4.85) 



(4.86) 
(4.87) 



Z£ 



Z* 



m 



z 6,c=l 

C+n { - n + (2m| o + m x + 2m| + mj)^ 

-\ m j a + 2m iS m x ~ m z + m2 t) + m %{ 2m x + m Z + ^iWl 
+{ml + m 2 i ){m} a -m 2 z ) 2 ^] 

+C? U • 2m e m x {T 2 t - 2(mj a + m|)7? + (mf^ - m|) 2 T '} 



(4.88) 
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neutralino (xf) exchange: 



w 



(x°) 
zi 



i 4 
mi . z — ' 



A °A 



»J=1 



C-,{T 3 « + ( a - m| - 2m 2 - m|)T 2 « 



+ [-s(m 2 + 2m|) + m x + m^(2m~ + 2m| - m 2 ) + 2m 2 m 2 z \T{ 



* y la 

-(mf o - mf^m 2 - m|)(m 2 + 2m|)T " 
+m x (m x o + m x o) ■ 3m 2 z [T^ - (m~ - mj)Tg 



+m x om x o [(— s + 2m z + m e )T^ 
+{s(m 2 + 2m|) - m|(m| + m 2 + 2m|) - m| (m 2 + m|)}T n ]| 
+C^-i • m^6m x m|77 + (m x o + m^p? - (2m 2 - m 2 z )T? 

+{m 2 x - m 2 z )(m 2 x + 2m|)7^] + 6m x m x om x om|7^| 



• lepton (£) - slepton (4) interference: 



w 



(i-i) 
zi 



1 A C^ z 
— ^-Re > 



m 



7 — s — m 
z 6=1 * 



+ii; j i 



where 



J 2 



— s 2 — s(m~ + m 2 + mfz — 2m~ — m 2 ) + (m~ — m 2 )(3m 2 z — m 2 ) 
+2|s[m~ — m~ (m| + m 2 + m|) + m 2 (m| — m 2 z )\ 

, o 2 r 2 / 2 2 \ , 2 / 2 , 2 2 \i 

+2m z [mj [mj - m x ) + m x (m x + m z - mj )J 
-m 2 [(m| - m~ a )(m~ a - m 2 x ) + m|(m~ + m 2 )]}.^*, 
s + 2ml + m 2 — 3m| — mfz — m 2 (m?7 — m 2 )(m~ — m 2 ) 

lb la g la ^ 

+2 js(m| - m| o + m|) + (m| - mfj 2 - m|(m| + mf) 
-m~ (m| - mf) - m 2 z (2m 2 x - m 2 )}.? 7 *, 
li* = -s - m 2 z + mj + - m| - m|) - (m| - m 2 z )(m 2 z - mf) 

—m~ (m 2 z + m 2 ) + 2m|m 2 j.7 7 *; 

lepton (£) - neutralino (%°) interference: 



w 



(i-x°) 

Zl 



i=i 



s — 



(4.89) 



(4.90) 



(4.91) 



(4.92) 
(4.93) 
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££Z*s-iaa 



iUZ* r^aa 



rZi 



fsyMZ*s~taa ^iUZ* r\aa i /~iMZ* s-iaa s~tUZ* r^aa \ r» , jZ 



iUZ* r\aa 



,2 jZl 



UZ*^iaa 



~i££Z* r\aa 



rZi 



where 



T Zi 



s 2 + s(2m\ + 3m| — m~ — m 2 — m 2 ) — (m~ — m 2 y ){3m 2 z — m 2 ) 
+2[s 2 m x o(m x + m^) + S [-2m|(m| - m 2 ) 



+ (m^g - m^)(m^ - m~ ) — mj(m o + m x ) 2 + m 2 z m o(m x o + 



-ra~ [2m|m?, + m?(m^ + 2m|)] 



J 5 



-to 2 (2to| - mj)(m~ - to 2 ) + to|to 2 (to 2 + m 2 z + m 2 ]^T u , 
m x o(—s + 2to| + to|) + 3m 2 z m x + |.sto x o(to 2 + 2to| 
+m x o [to|(2to 2 o - to 2 - to 2 - 2to|) + to 2 (to 2 - 
+3to x to|(to 2 — m~ + m 2 )\T u , 

Aj -t a J 



J 7 



s - 2to 2 - m~ + 3m?, 

A-i £a a 



3m 2 z 



-(m| - mj)(mj 



— 2{(m x o + to x ) 2 — m 2 z }{(m x o — m x ) 2 + 2m 2 z )T u , 
m x (s + m 2 z — m 2 ) + js[m x (TO 2 — m 2 + m 2 z ) — 3m 2 z m x o] 
+m x [m 2 (m 2 z - to 2 ) + m 2 z (m~ - 2m 2 + 2m 2 z ) + to 2 (to 2 - 2m 2 z )} 



+3m x om 2 z (mi a - m 2 )}^ 1 ; 
• slepton (£b) - neutralino (x®) interference: 



w 



Zl 



-x°) _ 



TO' 



7 .4 



? 6=1 i=l 

£>+ H {^3 - + m| o - to 2 + to| + mj 

n 2 x )]y 2 + ii"y 1 + i^y } 



+ -{m 2 z -mj)(mj - 

s Za 

-D%m e (m x o - m x ){ - 23^2 + (m| - to 2 )^i 
+ [s(m~ + to, 2 - 2to|) - to- + to| (3to| - to 2 ) 
-m 2 (to 2 - 3to|) - to, 2 (to, 2 + 2m 2 z )]y 

+ 1(TO| - TO 2 )(TO| - TO 2 )[^1 - (TO2 - ^ )3 ; o] 



(4.94) 



(4.95) 



(4.96) 
(4.97) 



(4.98) 
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._L (m |_ m 2 )2(m |_ m 2 )2;) ; o j 



where 



(4.99) 



/f = m x m x0 ( s + m | - m *) + >i (2m 2 - 3m|) 



2 ^ 



1 2 2 , 1 2/ 2 2 , r, 2 \ 

--m z m x + -m e {m~ a -m x + 2m z ) 



+ Ys ijnz ~ m ^ m I a ~ m x^ m I a -™ x + m z + m i) 

+4^ (m ^ ~ m ') 2 K a - m x) 2 > ( 4 - 10 °) 

-s[m~ (3m| + 3mf — 2m 2 ) + m x (m z + mf) — 6m$m z ] 

-m~ (m 2 x + m| + 2mf ) + m| {m,^ + m|m 2 - m z + mf (ra 2 + 4m| - 2m 2 )} 

—m z {m x + m z m x — m|(4m^ — m| — mf)} 



— — (m| — mf)(m~ — m 2 ){2m~ + m~ (8m 2 — 3m z + 7mf) 

+2m^ + 3m z m x + 2m| + (m 2 + 8m| + 2m|)} 

— m x m x o(s + m~ — m 2 ){(s — mf) 2 — 2m|(s + mf) + m^} 



2\2/ 2 2 1 2 1 2\ 

m ) (m~ -m Y + m z + m £ J 



2\3 



(4.101) 



4.22 £ aX -> 7^ 

This process proceeds via the s-channel lepton (£) exchange, and the t-channel slepton 
(£ a , a = 1,2) exchange 



~(£) . ~(e) . 
w t aX - n t = w fi + U V +w it 



lepton (£) exchange: 



w 



-it 



Cl a u (s 2 - 6mfs + mf) jl - ^(mf - m 2 ) 



(s — mf) 2 
-Cfn • 4m e m x (s + mf ) 



(4.102) 



(4.103) 



slepton (£ ) exchange: 



= 2e 2 



C^ u {7? + (m| - m 2 - mf)T/ - (m 2 + mf)™ 2 ^} 



-C-n^m^lrZ + mf^} 



(4.104) 
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lepton (£) - slepton (£ a ) interference: 



w 



(e-e) 



2e 2 



s — m 



r Re 



Cl a u js + 2m 2 - 2mf - 2mj 



+[s(mj a +m 2 x + mj) - 2(m£ - m 2 ) 2 - mj(m^ + m\ - mj)}^} 
+C™ n ■2m e m x { 1 + (« + 2m| -2m 2 . (4.105) 

4.23 £ aX -> W~v t 

This process proceeds via the s-channel lepton (I) exchange, the ^-channel sneutrino (uf>) 
exchange, and the it-channel chargino (xt, k = 1,2) exchange 



~(jt) . ~(V) , -(x*) , , -{l-x*) , -(^-x*) 



lepton (£) exchange: 



w 



nrj(.V(W+\2 



mly(s - m 2 ) 2 



(s - m 2 v )(s + 2m w ) 



2 m 2 
H 1 



sneutrino (ug) exchange: 



w 



(V) 



m 



w 



V*,taW+ (jxvm 



T 3 * + (m 2 + 2mj a + 2m 2 w )Tj 



{(ml - m 2 w ) 2 + 2m\(m 2 + m^)}7? + m 2 (mf_ - m^) 2 ^ 



chargino (x^ ) exchange: 



(4.106) 



(4.107) 



(4.108) 



(x±) 



£ c xtl a ve* c xf^ C ^ xW '*C^ xW '{t^ + (s- m\ a - 2m\ - m 2 w )T 2 u 



W k ,l=l 

+ [-s(m 2 + 2m 2 w ) + m 2 (2m| + m\ + 2m 2 v )]T 1 u 



-m~ l (m x - m w )(m x + 2m w 
+C^ W -*Cp W -m xt m xf { - (s - 2m 2 w )T? 



+ [s{m 2 x + 2m 2 w ) - m 2 w (m\ a + m 2 + 2m 2 v )]T M } 

+ m x (m xt C^ W -* C f* W - + rn xf C^ W -*Cp W -) • 3^{ - T« + m| 



(4.109) 
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lepton (£) - sneutrino (i>g) interference: 



w 



(t-u) 



2 (jlv ( w + *(ji>* l t a w + (jxyivi 
-^Re s £ 



m 



s — 



s{m~ - m~(m~ + m 2 + m 2 ^) + mj(m~ - m 2 ^)} 



+2m^{m| (m~ - m 2 ) + ra 2 (m 2 + - m-jlj.? 7 * j 
■ 2m e m x {s + + s(rn% t + ntfy - m%J 

+m 2 v (m 2 ^ + m~ - 2m 2 - m^)J JP"*| 



lepton (£) - chargino ) interference: 



(4.110) 



-x±) _ 



2 2 

— H,, X! - — 



s — m« 



+2 



c xe a e c xU w *f s 2 + s ( 2m 2 ± +Sm 2 v - mj - ml) - 3m^(m? - m 2 ) 
s 2 m 2 ± + s{(m 2 ± - m 2 )(m 2 ± - m~ ) + m 2 v (m 2 + - 2m~ )} 



+2m^(m| - m 2 )(m| - m 2 ±)]^"} 
+C x« c x^-*_ 2m ^| s _ 2m ^ 

+ s(m 2 ± - m 2 - 2m 2 w ) + m^mf^ +m 2 x + 2m 2 w - 2m 2 ±)].T t } 

_ c fai c xtxw-* 2mxm ^( s _ m ^ )(s + 2^)^ 
sneutrino (i^) - chargino (xt) interference: 



(4.111) 



w 



(?-x±) 



-Re p" ( *<«M' + *^x^''^^^« l/ ' 



m 



s + m~ - m^+rnly + -m 2 v (m 2 - m 2 ) 



3>2 



l r 

+ 2 



m 2 t (2m 2 x - Zmly) - rn^rn^ + ^m\ v {ra 2 ^ - to 2 ) 2 



+ -m 2 w {m\ a - m 2 x ){m\ a - m 2 + m 2 w )\ y 1 + J^o} 
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+Cft k xW *m x m x ± | - (s + m^)3^i + (s - mwY 



l + -(mj -m 2 ] 



y } 

(4.112) 



where 

jWu 



-s[m~ (3m^ - 1m\) + m 2 w m x \ 

4/2, 2 \ , 2/4, 2 2 4 \ 2 2/2, 2 \ 

-mj (m x + m w ) + m~^(m x + m w m x - m w ) - m w m x (m x + m w ) 

~^ m w( m q a ~ m x)^ 2m f a + m f a ( 8m x ~~ 3m w) + 2m x + 3m^m^ + 2m^/} 
__ m ^( m | _ m x ) 2 (m| - m x + m^) - — m^(m| - m x ) 3 . 



(4.113) 



4.24 -> 



The process £ a x — ► ^ involves the s-channel lepton (£) exchange, the i-channel slepton 
(£ a , a = 1,2) exchange, and the u-channel neutralino (%°, i = 1,2,3,4) exchange 



™J aX ^« = w ht + w hi + W M + W M + W M + W M 



(4.114) 



lepton (£) exchange: 



w 



he 





2 " 


s — m 2 





4Ct a n mem x (2s + 2mj - ml) 



+C™ u {s 2 - s(m 2 h - 6m 2 ) - m 2 {m\ - mf)}{l - -{v^ - m 2 )} 



(4.115) 



slepton (£;,) exchange: 



4} = E Cf ah cf° h * \CX\A-Tl + K + mfrt} + 2C b J n mem x ^ 



b,c=l 



(4.116) 



neutralino (xf) exchange: 



w 



(\") \- r vV Y A ' v " 

/ ^ S S 



Id 



C%i{ -(s-m 2 x - mj)T- - (m| - m 2 )(m 2 x - m 2 )T u 



+m x (m x o + m^p? - (m| - m 2 e )T u ] + m^om^T? 1 + (a - mf^ - m|)T u ]} 

+C^ i m^|2m x T 1 u + (m x o + m^pj* + (m£ - m|)T u ] + rn^m^o • 2to x T m } 

(4.117) 
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lepton (£) - slepton (4) interference: 



w 



{t-t) 
hi 



2Re E 



s — m„ 



C? n m x (s - m£ + 3m|)JT' 



2\ -r-t 



+Cf u m e [ -l + (s + 2m 2 x + mj - m| - mf ; 
lepton (£) - neutralino interference: 



(4.118) 



hi 



4 r^Uh*r< x i xh 
2Re^ ° 5 ° 5 



i=i 



s — m; 



+ [2m x os + (m x o + m x )(m x o + m 2 x - 2m~ - ml) + 2m x m 2 ]J 7U ^ 
+C""j|s + m 2 + [sm^fm^o + m x ) — m x m x o{n? h — 3m 2 ) 
-mjml + mj(m 2 x o + 2m 2 x - ml)]?"} 

slepton (4) - neutralino (xi) interference: 



(4.119) 



2 4 



6=1 i=l 



+ 



— m x o(s — m\ + mf){l — 



1 

-{m 

s 



(mj - ml)} 



l | o o o o 

+C_ li me< yi + s — m~ + m ft — 3m x — m e — 4m x m x o 



+m x {s + m 2 l -m 2 h -m 2 x -2,m 2 --(m 2 h -m 2 )(m 2 l - m 2 x )}]y Q ) 

(4.120) 

The contribution from £ a \ — > can be obtained by a simple substitution h ^ H. 
4.25 £ a x -> A£ 

This process involves the s-channel lepton (£) exchange, the ^-channel slepton (£ a , a = 1,2) 
exchange, and the u-channel neutralino (%°, z = 1,2,3,4) exchange 



• lepton (£) exchange: 



1 



w 



Ai 



= W At + W V + ^Ju + W Ai + w Ai + W Ai ■ 



4Cl a n mim x m 2 A 



(4.121) 



qUA 


2 " 


s — m 2 
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+C™ n {{s - m 2 ) 2 - m\{s + m 2 )} {l - \{m\ a - m 2 x )} 



(4.122) 



slepton (£},) exchange: 



w 



(0 _ n t* b t a A e*e a A* 
ai — P P 

b,c=l 



CWA-Tl + {ml + mj)Tj} + 2C b _ c n m e m x ^ 



(4.123) 



neutralino (xt) exchange: 



(x°) 



cf iXA c x ) xM 



C%.{ -{s-ml- m 2 )^ - {ml - m 2 ){m 2 x - m 2 A )T u 
-m x {m x o + m x o)[Tf - (m| - mf)T M ] + m^m^T? + {s - m| - m 2 A )T^]} 

-C^me^m^ - {m x o + m x o)[Tf + {m\ - m 2 A )T£] + m^m^o • 2m x T u } 

(4.124) 

• lepton {£) - slepton (£(,) interference: 



w 



(t-t) 
Ai 



2Re E 



2 riUA*r i ^ a A 



8 -mi 



D a _ b u m x {s -m\- m 2 )P 



+D? 11 m e { - 1 - {s + m| - mf a - m £ 2 )^'}] ; 
• lepton {£) - neutralino (x°) interference: 



(4.125) 



w 



Ai 



2Re E 



Op Op 



s — 



C™ u me^m x o — m x 

+ [(m x o - m x ){m 2 x o - m x ) - {m x o + m x )m^]J rU | 

~~ s + m 1 + [ — sm x o{m x o — m x ) — m x m x o{m 2 A + m 2 ) 



+mjm 2 A + m?(m x o - m^-T 1 } 
• slepton (4) - neutralino {x®) interference: 



(4.126) 



2 4 



it; 



6=1 i=l 



- 32 - 



+ 



— m x o(s — m 2 A + m|){l (m~ — m^.)} 



-m x { S + mf n - m\ - m\ - 3m 2 - ±{m 2 A - m, 2 )(m| - m 2 )}]^} 

i I 

— D_ li mi<.yi + s — m~ + m A — 3m x — m { + 4m x m x o 

-i(,4- m ?)K.-^)]^} j 

(4.127) 

4.26 -> 

This process proceeds via the s-channel lepton (^) exchange, the t-channel sneutrino (z/^) 
exchange, and the tt-channel chargino (%^, fc = 1,2) exchange 

wj „_ =w { S- + w%L + +w^+w { "Z x±) +w^: x±) : (4.128) 

• lepton (£) exchange: 



w 



H-vi 



( jv l tH+ 


2 

(s 


s — m 2 




al 2 



(s - m 2 H± ] 
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chargino (Xk) exchange: 



(4.129) 
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lepton (£) - chargino (xt) interference: 
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(4.134) 



5. Numerical Analysis 

In this Section we present some numerical examples to illustrate the effect of the neutralino- 
slepton coannihilation. We will work in the framework of the CMSSM where the effect 
considered in this paper is particularly important. In computing the neutralino WIMP 
relic abundance we will employ the exact expressions for the cross-sections of neutralino 
pair-annihilation derived in jl and the neutralino-slepton coannihilation ones listed above. 
To generate mass spectra we will use the package SUSPECT (v. 2. 05) which includes 
full one-loop radiative corrections to sfermion masses as well as to the effective potential. 

We begin by presenting in Fig. |l| contours of the relic abundance f^/i 2 in the plane 
(mi/2, mo) for tan (3 = 10 (left window) and 40 (right window), and for Aq = and [i > 0. 
The solid (dashed) curves correspond to including (neglecting) the coannihilation effect. 
The light orange region of fl x h 2 > 0.3 is inconsistent with the age of the Universe while 
the green band corresponding to 0.1 < £\/i 2 < 0.2 is favored by direct measurements of 
the dark matter component in the Universe. For the sake of clarity, we only denote (in red) 
the regions of the plane where the lighter stau f\ is the LSP (and in some part of it one of 
the Higgs boson mass-square is negative), as well as those excluded by the LEP limit on 
the lightest chargino but no other experimental bounds. In particular, we do not indicate 
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Figure 1: Contours of £l x h 2 in the plane (m 1 / 2 , mo) for tan/3 = 10 (left window) and tan/3 = 40 
(right window), and for A = 0, /i > 0, m v t olc = 175 GeV and m fc (m 6 )|f = 4.20 GeV. The red 
regions bands are excluded by chargino searches at LEP and corresponds to the lighter stau being 
the LSP. The light orange regions of fl x h 2 > 0.3 are excluded by cosmology while the narrow green 
bands correspond to the expected range 0.1 < il x h 2 < 0.2. Also shown are the semi-oval contours 
of fl x h 2 in the absence of coannihilation. 



the regions inconsistent with the lightest Higgs boson mass bound from LEP, nor with 
BR(£> — > X s j). (Their effect has been presented in with updates in |$^].) We only note 
that these bounds, if taken at face value, exclude the favored regions of 0.1 < Q x h 2 < 0.2 
(green bands) of m 1/2 < 320 GeV (tan/? = 10) and m 1/2 < 550 GeV (tan/? = 40). In 
reality, theoretical uncertainties may considerably weaken these bounds, but in any case, 
it is clear that the remaining region is allowed mostly only by the coannihilation effect. 

Note that the narrow band opened up by coannihilation eventually ends at the bound- 
ary of equal neutralino and stau masses. This is in qualitative agreement with the results 
obtained in |l|, |23|, |24 ] using the usual partial wave expansion but not with the more recent 
analysis [16| where the cosmologically allowed region appears to lie basically parallel to the 
boundary even at very large m^. We will come back to discussing the upper limit on m x 
below. 

The effect of coannihilation becomes dramatic when the two coannihilating particles 
are nearly degenerate in mass, but can be significant even for the mass difference of some 
20 — 40 GeV. This is illustrated in Fig. ^ where the relic abundance £l x h 2 is plotted as a 
function of Am = m x — for tan (3 = 10 (left window) and 40 (right window) at a fixed 
value of vn-yji = 500 GeV. Also marked are the cosmologically excluded (f2 x /i 2 > 0.3) and 
favored (0.1 < & x h 2 < 0.2) regions. 

In Fig. H we present the total as well as the individual contributions to the quantity 
J(xf) = Jq 1 dx(avy[0\){x) from the various classes of processes listed in Table p]. (Roughly, 
Q x h 2 ~ 1/J(xf) - see, e.g., [f|.) This is done for a slice of constant mi/2 = 500 GeV as a 
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Figure 2: 



The relic abundance fl x h 2 



with (solid) and without (dash) coannihilation vs. Am 



m x — m^ 1 for tan/3 = 10 (left window) and 40 (right window) at a fixed value of mi/2 = 500 GeV. 
Also marked are the cosmologically excluded (il x h 2 > 0.3) and favored (0.1 < il x h 2 < 0.2) regions. 
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Figure 3: Total and partial contributions to J(xf) from the various classes of processes listed in 
Table MS iX function of Am/m x = (m x — mj l )/m x for the two choices of tan/3 = 10,40 of Fig. |l|. 

function of Am/m x = (m x — m^ 1 )/m x for the two choices of tan/3 = 10,40 of Fig. |l|. At 
a smaller mass difference it is the slepton-slepton annihilation into lepton pairs (mostly 
t\t\ — > tt) that appears to be dominant but it drops quickly and at larger Am/m x it 
is overtaken by the neutralino-slepton coannihilation channel (mostly the 7T final state). 
Here our conclusions qualitatively agree with pi. 



As mentioned in the Introduction, [23| contains the only set of analytic expresssions for 
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tan0=lO. A D =0, ,ii>0 



tan0 = 4O. A D =0, fj.>0 





500 600 

m y (GeV) 



500 600 

m, (GeV) 



Figure 4: The minimum of fl x h 2 (i.e., fl x h 2 along the line m x = m^) as a function ofm x for the 
two representative choices of Fig. |l| The solid (red) curve correspods to the exact neutralino-slepton 
coannihilation cross sections, while the dashed (blue) ones to the approximate ones of P3|, 



the neutralino-slepton coannihilation that is at present available in the literature. As noted 
by the authors, the expressions given there are approximate as they were derived using 



partial wave expansion. They also become less reliable at large values of tan/3 > 20 [24] 
because the effects of h T and of the t\ — 75 mixing, and in some channels of the mass of 
the r, were neglected. We have made an attempt at improving the expressions of [^] by 
including in the propagators the widths of the gauge and Higgs bosons and the neutralinos; 
otherwise they become singular. 

While keeping these points in mind, it is nevertheless interesting to compare them with 
the exact expressions derived here. We do this numerically in Fig || for the cases displayed 
in Fig. [j]. Along the boundary of m x = (and in fact in most of the (mi/2i m o) plane) 
m x ~ 0. 44771^2 — 2.8 sin 2(3 |33| j. This allows us to plot the minimum value of £l x h 2 along 
the boundary as a function of m x . It is clear that £l x h 2 increases with m x and at some point 
becomes inconsistent with Q x h 2 < 0.3. The solid (red) curve correspond to the exact results 



while the dash (blue) ones were obtained by using the approximate expressions of [23]. For 
tan/3 = 10 the upper limit from Q x h 2 < 0.3 on m x is strengthened from ~ 700 GeV down 
to some 640 GeV. In contrast, at large tan/3 = 40 (where, we repeat, the approximate 



expressions are not really applicable [23]), an upper bound becomes considerably weaker 



(m x < 765 GeV) than what one would get by naively applying the expressions of [23]] . The 
ranges of m x corresponding to the favored range 0.1 < £l x h? < 0.2 also shift accordingly. 



6. Summary 

The accuracy of determining the abundance of the dark matter in the Universe is continu- 
ously improving. This requires theoretical computations of the neutralino relic abundance 
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to be performed with at least the same, if not better, level of precision, if one wants to 
reliably compare theoretical predictions with observations. 

In this paper we have derived a full set of exact, analytic expressions for the neutralino- 
slepton coannihilation cross sections into all tree-level two-body final states. While these 
formulae are applicable in the framework of the general MSSM, they are of particular 
importance in the context of the Constrained MSSM. In this framework, which is often 
considered a "reference" SUSY model and is thus of much interest to the community, much 
of the allowed regions are a result of the neutralino-slepton coannihilation. Our results 
should help in allowing one to determine these regions more precisely. 
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A. Lagrangian Terms and Couplings 



In this Appendix, we define the couplings which appear in the main text and which have 
not been defined in Hj. We follow the same conventions and notation as in that paper. 
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Some more couplings including the neutrinos: 
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The couplings £ x i iL , £ x i iR , r x ^ L and r x ^ R are defined in 0]. Note that I = e, //, r 
represents a charged lepton. We neglect generation mixing in the lepton sector. The 
charged slepton mass eigenstates i a (a = 1,2) are related to the slepton gauge eigenstates 
£l and £r via 

la = (V t ) al £ L + (VtWR, (A.46) 

where Vi denotes a 2 x 2 matrix which diagonalizes the charged slepton mass matrix: Vi 
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B. Auxiliary functions 



Here we give definitions for the auxiliary functions used in the text. The functions T 1 and 
T u are given by 
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For yi = y 2 , D\ vanishes so that the two functions and F u reduce to the same function 
which is denoted by T . 

The functions and 7^" are obtained from % in Q by the following replacements 
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The functions % in the right-hand side represent those defined in Q , where F in Q should 
be replaced with F in this appendix. For y x = y 2 , the two functions and 7^" reduce to the 
same function which is denoted by %. The expressions for the functions 3^ (i = 0, 1, 2, 3, 4) 
are given by 
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&c 




£ J* _> y A 

0> h ! 










7 7* _> hh f{H hH 


h,H 


-^c 




4P 


7 7* — > h A HA 


A,Z 








7 7* _> 4 4 


h,H 


■^c 




4P 


77*^ h+ ff- 


h,H,-y,Z 




vt 


4P 


7* i />i 


h,H,A,-y,Z 


Xi 






— > a a 

^ o ^ 5 


h,H,A,-y,Z 








4^ -> ^ 


Z 


x^ 






7 7'*^ fp 




x^ 






4ift* W 














x° 






44 ^ ^ 




x° 






<-OA Zyt 


e 




x° 




7 Y _> /-y/ 


e 








4x - 


e 


^£ 






IaX -► 


e 


&c 


Xi 




4x - ^ 


e 




x° 




4x 


e 


7e 


Xfc 





Table 1: A complete set of processes relevant for the neutralino-slepton coannihilation into tree- 
level two-body final states in the MSSM. 'PI' denotes four-point (4P) interactions. The notation 
is as follows: £ a denotes e a , Jl a and r a , where a = 1,2 (and likewise b and c) denotes a slepton 
mass state index for each generation. The symbols I and £' represent charged leptons of different 
generations, and likewise for the sleptons. For the neutralino index i = 1, . . . , 4 and for the chargino 
one k — 1,2. Note that for each reaction with a non-zero net electric charge there is a corresponding 
one with the the opposite net charge which is not included in the Table. For example, in addition 
to Ij b -> it there is also I*J* b -> 11. 
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